Hidden nonlinear supersymmetry of finite-gap Lame equation 
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Supersymmetry as a fundamental symmetry pro- 
viding a natural mechanism for unification of gravity with 
electromagnetic, strong and weak interactions, still waits 
for experimental confirmation. On the other hand, in 
nuclear physics supersymmetry was predicted theoreti- 
cally 0] and has been confirmed experimentally [sll as 
a dynamic symmetry linking properties of some bosonic 
and fermionic nuclei. It would be interesting to look for 
some physical systems whose special properties could be 
explained by a hidden ordinary (not a dynamic) super- 
symmetry. 

In the present Letter we show that the quantum sys- 
tem described by the finite-gap Lame equation possesses 
a hidden supersymmetry. A very unusual nature of the 
revealed supersymmetry is that it manifests as a non- 
linear symmetry of a bosonic system without fermion 
(spin) degrees of freedom. This means that we find here 
a kind of a bosonized supersymmetry giving a natural ex- 
planation for peculiar properties of the periodic quantum 
problem underlying many physical systems. 

The Lame equation first arose in solution of the 
Laplace equation by separation of variables in ellipsoidal 
coordinates and one of its early applications was in 
the quantum Euler top problem [5[. It plays nowadays a 
prominent role in physics appearing in such diverse the- 
ories as crystals models in solid state physics [1, [v L ex- 
actly and quasi-exactly solvable quantum systems a , 
integrable systems and solitons [1^, su per symmet- 
ric quantum mechanics 12l. BPS monopoles [13| , instan- 
tons and sphalerons [IJ, |l5| , classical Ginzburg-Landau 
field theor y [l^ . Josephson junctions [o], rnagnetostatic 
problems [18|], inhomogeneous cosmologies [19|, Kaluza- 
Klein theories [13] , chaos 21 1 , and preheating after infla- 
tion modern theories p^ . Most often, the Lame equa- 
tion appears in physics literature in the Jacobian form 
of a one-dimensional Schrodinger equation with a doubly 
periodic potential, 

H,^ = E^, H,=-^+j{j + l)k^sn^x,k), (1) 

where sn(x, k) = snx is the Jacobi elliptic odd function 
with modulus fc (0 < A; < 1), and real and imaginary 



periods AK and 2iK', K — K(k) is a complete elliptic 
integral of the first kind, and K' = K{k'), fc'^ = 1 - fc^ 
0, ■ ^ remarkable property of this equation is that at 
integer values of the parameter j = n, its energy spec- 
trum has exactly n gaps, which separate the n-\-l allowed 
energy bands. The 2n -t- 1 eigenf unctions associated to 
the boundaries Ei{n), i — 0,1, ... , 2n, of the allowed en- 
ergy bands [Eq, Ei], [£^2, £^3], ■ • ■ , [E2,i, 00] are given by 
polynomials ('Lame polynomials') of degree n in the el- 
liptic functions snx, cnx and dux. These polynomials 
have real periods AK or 2K, and the boundary energy 
levels Ei{n) are non- degenerated. The states in the inte- 
rior of allowed zones are described by the quasi-periodic 
Bloch-Floquet wave functions (which can be expressed in 
terms of theta functions ^) of quasi-momentum k{E), 

*|(a; + 2K) = exp{±iK{E))^^{x). 

Every such interior energy level is doubly degenerated. 
For any non- integer value of the parameter j, Eq. ([1]) 
has an infinite number of allowed and prohibited zones. 

The double degeneration of the energy levels is typical 
for a quantum mechanical system with N = 2 supersym- 
metry. But the presence of 2n -I- 1 edge-bands singlet 
states in the n-gap Lame equation indicates on an un- 
usual, nonlinear character [2^ of a possible hidden 
supersymmetry. To reveal it, one notes that in the lim- 
iting case k — 1 we have K — 00, K' — ^ , and system 
([T|) reduces to the Poschl- Teller quantum system given 
by the potential 

U{x) = -I- 1) sech^a; + j{j + 1). 

The latter, as it was shown recently in [2^, at j = n 
possesses a hidden polynomial supersymmetry 24, [2^ 
of order 2n -\- 1 generated by the supercharges Qn and 

[Qn, Hn] - [Qn,Hn] = 0, {Qn, Qn} = 0, (2) 
Qn = Qn = P2n+l{Hn), (3) 

where P2n-i-i(-ffn) is some polynomial of the degree 2n-\- 
1 of the Hamiltonian Hn, i? is a reflection, i?^'(a:) = 
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^{—x), identified as the grading operator, 

[R,H„]^0, {R,Qn} = {R,Qn} = 0, = (4) 

and Qn is a self-conjugate local differential operator of 
degree 2n + 1. Based on this observation, first we note 
that in the trivial case of the free particle system with 
j = characterized by the one allowed ('conduction') 
band [Eo{0), oo], Eo{0) = 0, the odd first order differen- 
tial operator Qq = —iD, D = is identified as the su- 
percharge. For the one-gap Lame system ([1]) with j ~ I, 
let us look for the self-conjugate integral of motion Qi, 
[Qi,Hi] = 0, in the form of the third order differential 
operator. A direct check shows that 

iQi^D^ + fD+^f, (5) 

is the odd integral, {R, Qi} ~ 0. Here 

f ■.= l + k^ -Sk^an^x, (6) 

/' — Df. The double-periodic elliptic function / with 
periods 2K and 2iK' satisfies the elliptic curve equation 

= IK -/)(/-a2)(/-a3), (7) 

whose characteristic roots are 

ai = /(0) = l + fc2, a2 = /W = l-2fc2, 

as^ f{K + iK')^k^ -2, (8) 

ai + a2 + a3 = 0. Differentiation of Eq. ^ gives the 
identities 

f" + 2f = 2b^ D\D^ + 2f)f = Q, (9) 

where &^ = —^(0102 + aioa + a2a3) = fc* — fc^ + 1, 
I = 1,2,.... Using these relations, one finds that = 

PsiHi) - {Hi - Eo{l)){Hi - Ei{l)){Hi - ^2(1)). 
The energies 

£;o(l) = A:2, £;i(l) = l, S2(l) = l + fc' (10) 

correspond here to the eigenfunctions ^E'q^'' dnx, 
^'^^^ = cnx, 5*2^^ = snx, which form a zero-mode sub- 
space (kernel) of the supercharge Qi. The states in the 
interior of the two allowed zones are described by the 
quasi-periodic eigenfunctions 

*1 = Q^^^ exp{TxZ{a)), (11) 

where H{x), Q{x) and Z{x) are the Jacobi Eta, Theta 
and Zeta functions, while the parameter a is related to 
the energy eigenvalue E via the equation E = dn^a + 



k^, see Ref. They are also the eigenstates of the 

supercharge, 

QiVl/± = ±v/7^vI/± (12) 
Assuming that the j = 2 Lame polynomials ^'J, = / + 

(2) f2) f2") 

b, = cna;dna;, = snxdnx, ^'3 = sna;cnx, 
^4^^ = / — ^ form a zero-mode subspace of the fifth 
order integral Q2, one finds 

iQ2 ^D^ + 5fD^ + ^fD^ + (§/" + 4/2) D. (13) 

In the same way for j — 3 and j = 4 a tedious cal- 
culation gives the supercharges Q3 and Q4. We do not 
display their explicit form here, but, instead, describe 
the general structure of the supercharges corresponding 
to j = 0,1,2,3,4. First, one notes that if the deriva- 
tive is assigned the homogeneity degree dh{D) = 1, in 
accordance with Eq. ^ the function / can be assigned 
dhif) = 2, and then dh{Hj) = 2 and dh{Qj) = 2j + 1. 
Every supercharge has the leading term D^^~^^, the next 
term is of the form fD^^~^, and every subsequent term 
decreases the order of the derivative on the right in one 
unit. The supercharges corresponding to the even cases 
j = 0, 2, 4 contain the last term to be proportional to D. 
With this structure of the supercharges for the first cases 
of J = 0, . . . , 4 at hands, we can fix now the form of the 
supercharges in the generic case j = n. Let us present 
the Hamiltonian operator in terms of the function /, 

H, = -D'-h,{f{x)-f{0)), h,^^j{j + l), (14) 

and look for the supercharge in the form 

+(5j/2)D2,-3 ^ A,/"'i?2j-4 + ...^ (15) 

where in coefficients associated to the factors I > 0, 
it is necessary to include all the independent structures 
of homogeneity degree dh = 2j + 1 — I given in terms 
of / and its derivatives modulo identities ([9]). Requiring 
[Qj, Hj] = 0, we can fix the first coefficients, 

aj = hj {j + i) , I3j = aj {j - i) , 7j = f/3j(j - 1), 
= U - 6) , A, ^-^^^{j-D ij - 2) . (16) 

These coefficients allow us to find a recurrence relation 
for the supercharges. We note that the kernel K2n of the 
supercharge Qj with j = 2n is spanned by the 4n -I- 1 
functions 

^„.(1, /,..., /"-I), a = l,...,4; r, (17) 

ipi = snxcna;, ip2 = snxdnx, 933 — cna;dna;, 954 — 1, 
which are linear combinations of the Lame polynomials 
of the degree 2n. For j = 2n + 1, the kernel 1^^2,1+1 is 
formed by the 4n -f 3 functions 

Xa-(l,/,...,r-'),a=l,...,4; Xar,a= 1,2,3, (18) 
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with xi = dnx, X2 = cna;, Xi = snx, X4 = 
smcnxdnx, where for n = the states proportional 
to X4 are absent. In comparison with the kernel Kj_2 of 
the supercharge Qj_2, the kernel Kj of the supercharge 
Qj includes four additional states. These are (paf^"^, 
a = 1,2,3, and /" for j = 2n, and Xaf"'^^, a — 1,2,3, 
and X4/"^^ for j = 2n + 1. Therefore, there should exist 
a relation 

Q, = AjQj_2, (19) 

where Aj is a fourth order differential operator of homo- 
geneity degree dh{A) = 4 of the form (modulo the first 
identity from 

A, =D^ + a^fD^ + p,f'D + 7,fe2 + f^f"_ (20) 

Using Eqs. p5|) . p6)) . one finds the numerical coefficients 
of the operator Aj 

a, - 2j{j - 1) + 1 , /3, = - 7) + I , 

7, = J-^(J- - If , - i(j + m + l)(j - If. (21) 

The operator Aj is not symmetric, and in the representa- 
tion it serves to annihilate the additional zero modes 
of Qj after application to them of the operator Qj-2- 
There is also an alternative recurrence representation of 
the supercharge, Qj = (5j_2Aj. The Hermitian conju- 
gate operator Aj acts invariantly on the kernel of Qj-2, 
A^ : Kj-2 and transforms four additional states 

of the kernel of Qj into some linear combinations of the 
states of Kj-2- Relation ^9]) (or, Qj — Qj-2^j) allows 
ones to calculate Qj for arbitrary even and odd values of 
j proceeding from the explicitly displayed supercharges 
Q2 (or, Qo) and Qi. 

With the fixed form of the integral Qj, let us discuss 
the general structure of a hidden polynomial supersym- 
metry. Since Qj is a self-conjugate odd local differen- 
tial operator, one can introduce another, nonlocal su- 
percharge, Qj — iRQj. Hj and R are commuting self- 
conjugate operators. Let be their common eigen- 
states, Hj^^ = R^e = Since Qj com- 

mutes with Hj and anticommutes with R, there exist 
some linearly independent combinations ^'_E,q and R"^ E,q 
of ^E"^ and ^'^ such that 

Qj-^E,q^q{E)-^E.q, QjR'fE.q^ -q{E)R^E,q- (22) 

Then, the states "^E.q and R^E,q are the eigenstates of 
Qj and Qj with the same eigenvalue q^{E), and, hence, 
the same is vahd for the states ^ff, Qj*^ = Qj^^ = 
q^{E)'^^. All the states ^'^ corresponding to the al- 
lowed zones constitute the basis in the class of Bloch 
functions (including the periodic and antiperiodic ones). 
Hence, we get the operator equality Q^ = Q^ = q^{Hj). 
The operator Qj is a local differential operator of de- 
gree 4j -f 2. This means that the operator q^{Hj) 



is a polynomial of degree 2j -I- 1 of its argument, i.e. 
q^{Hj) = C{Hj — co){Hj — ci) . . . {Hj - C2j), where C is 
a real constant, while Ci are real, or some pairs of them 
could be mutually conjugate complex numbers. Compar- 
ing the coefficients in Q j and {Hj ) before the operator 
we find that C — 1. Then, applying the opera- 
tor q^{Hj) to the 2j + 1 (anti)periodic eigenstates of the 
Hamiltonian Hj corresponding to the boundaries of al- 
lowed zones Ei{n), and remembering that the same states 
constitute the kernel Kj of the supercharge Qj , we find 
that the set of the constants Ci coincides with the set of 
the boundary eigenvalues Ei{n), i — 0, . . . , 2n. Thus, we 
have shown that Q^ = Q^ = P2n+i{Hn), where 

P2n+i{E):=UT:oiE-Edn)) (23) 

is the Lame spectral polynomial. The nontrivial odd in- 
tegrals generate the order 2n+l polynomial superalgebra 
being the hidden symmetry of the bosonic system ([!]). 

The nonlinear character of the local supercharges and 
supersymmetry of the system ([I} is reminiscent to a non- 
linear symmetry of a particle in a Coulomb potential gen- 
erated by the Laplace-Runge-Lenz vector integral, and 
to that of an anisotropic oscillator with commensurable 
frequencies [27| . Let us clarify the dynamical picture un- 
derlying the hidden nonlinear supersymmetry having in 
mind the analogy with the anisotropic oscillator. Con- 
sider the one-gap case. The Hamiltonian Hi can be fac- 
torized in three possible ways: 

Hi = A^Ad + P= AlA, + 1 = AlA, + l + k\ (24) 

where Ad = D — (In dnx)', A^dnx = 0, and Ac and As 
have a similar structure in terms of the cnx and snx. 
Write the Heisenberg equations of motion of Ad and Al, 

iAd^ LUd{x)Ad, iAl = -AIlUs{x), (25) 

ujd{x) — — 2(ln dna;)", uJs{x) = — 2(ln snx)". Define the 
operator Ag/d ~ D — (In snx)' + (In dux)' , for which 

iAs/d = ujs{x)As/d - As/dUJd{x). (26) 
Then the relation 

tAlAs^dAd = Qi (27) 

gives us one of the six possible factorizations of the super- 
charge dS]). Note that operators Ac, Ag, A^/d and associ- 
ated instant frequencies have singularities on a real line, 
which cancel in the Hi and Qi. For j = n > 1 the same 
dynamical mechanism underlies the supercharge struc- 
ture and its possible factorizations. 

We conclude that the physical systems associated with 
the n-gap Lame equation possess a hidden bosonized non- 
linear supersymmetry. It is behind the double degener- 
ation of the energy levels in the interior of the allowed 
bands and the singlet character of the 2n -I- 1 edge-bands 
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states. The latter form a zero-mode subspace of the local 
supercharge Qj (as well as of the nonlocal one, Qj) be- 
ing a differential operator of degree 2n+ 1. Taking into 
account the parity of the states p7)) and p^ . one finds 
that the system (IT|) with any j — n is characterized by 
the Witten index [28| equal to 1. The information on the 
transfer matrix can also be extracted from the structure 
of its hidden supersymmetry. The detailed analysis of 
this aspect will be presented in a separate publication. 

In the limit A; = 1, sna; = tanh x, ca.x — dna; = 
secha;, the valence bands [Eq^Ei], . . . ,[E2n-2,E2n-i\ 
shrink, and two boundary states of a valence band trans- 
form into one bound state of the related Poschl- Teller 
system. As a result, the kernel of the supercharges of 
the latter system is constituted not only by the bound 
eigenstates and the lowest eigenstate from the continu- 
ous part of its spectrum, but also should include some 
n unbounded states. The discussion of this limit of the 
Lame equation corresponding to the Poschl- Teller prob- 
lem, and its relation to the bound state Aharonov-Bohm 
and the Dirac delta potential systems from the viewpoint 
of the hidden supersymmetry will be presented elsewhere. 

A simple shift of the argument in the Lame equation 
for a half of the real period of the Hamiltonian , a; — > 
X + K, gives the isospectral doubly-periodic system 



= -D^+j{j + 1)(1 - dn-2(x, k)). 



(28) 



At j = n it has a hidden polynomial supersymmetry 
generated by the supercharges Qj and Qj , whose explicit 
structure can be obtained by applying the Jacobi func- 
tions identities sti(x + K) = cna;/dna;, cn{x + K) — 
— fc' snx/ dna;, dn {x + K) = k' / dux to the supercharges 
of the system ([T|) . It would be interesting to look for the 
hidden polynomial supersymmetry in other finite-zone 
double periodic quantum systems. 

Finally, it would be interesting to clarify the role played 
by the revealed hidden nonlinear supersymmetry of the 
finite-gap Lame equation in the periodic Korteweg-de 
Vrics equation theory |lll| and in the periodic relativistic 
field theories 14 , 1^ 16| to which system ([1]) is intimately 
related. 
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